x,x0,x1=[0,0.01,2*pi];
L1(x)=(x-x1)./(x0-x1).*y0+(x-x0)./(x1-x0).*y1;
y0=sin(x0);
y1=sin(x1);
plot(x,x0,x1,y0,y1,L1(x))
